30th International Cosmic Ray Conference 




ICRC*/ 



Merida, Mexico 



o 
o 

(N 
Oh 

00 



X-Ray Observations of LSI+61303 with Swift 

J. Holder 1 , A. Falcone 2 and D. Morris 2 . 

1 Bartol Research Institute, Department of Physics and Astronomy, University of Delaware, Newark, DE 
19716, USA 

2 Department of Astronomy and Astrophysics, Pennsylvania State University, University Park, PA 16802, 
USA 

jholder@physics. udel. edu 



43 
Oh 

6 



(N 
> 

00 

m 

(N 

ON 

O 

o 



X 



Abstract: The TeV emitting high-mass X-ray binary system LSI+61303 was observed with the Swift 
satellite from early September 2006 to early January 2007. Many of these observations were contempora- 
neous with TeV observations. The data consist of observations on 24 separate days with durations ranging 
between 700s and 4700s, and partially cover 4.5 orbital periods of the binary system. We present here an 
analysis of the 0.2 to lOkeV X-ray data from the Swift-XRT instrument. Contemporaneous optical data 
from UVOT are also available. 



Introduction 

LS I +61 303 is a high mass X-ray binary system 
at a distance of ~ 2 kpc, composed of a B0 Ve star 
with a circumstellar disc and a compact object. The 
observed radio through optical emission is modu- 
lated with a period (P = 26.4960 ± 0.0028 days) 
believed to be associated with the orbital period of 
the binary system [6, 2]. Periastron takes place at 
phase 0.23 and the eccentricity is 0.72 ± 0.15 (al- 
though see [7] for a recent re-evaluation of the or- 
bital elements). 

The detection of an extended and apparently pre- 
cessing jet-like radio-emitting structure led to the 
possible identification of LS 1+61 303 as a micro- 
quasar, with the emission generated through par- 
ticle acceleration in a relativistic accretion-driven 
jet [10]. More detailed AU-scale observations [4] 
strongly support a counter-model, where the radio 
emission arises from particles shock-accelerated in 
the interaction of a pulsar wind with the circum- 
stellar material. The absence of any features in 
the X-ray spectrum (e.g. thermal components due 
to an accretion disc) also support the pulsar wind 
model. However, Romero et al. [11] argue that 
the relative wind strengths are such that it is not 
possible to produce the simple elongated shape ob- 
served in the VLBA images, and that the gamma- 



ray lightcurve is more easily explained as the result 
of variable accretion onto a compact object. 

At higher energies, LS 1 +61 303 was associated 
with the COS-B source 2CG 135+01 [13], and the 
EGRET source 2EG J0241+6119 [8]. More re- 
cently, LS 1 +61 303 has been discovered to be a 
strong source (peak flux ~ 10% of the steady flux 
from the Crab Nebula) of very high energy gamma- 
rays (> 100 GeV)[l, 9]; one of only three known 
galactic sources of variable TeV emission. The 
particle acceleration and photon generation mecha- 
nisms which produce the TeV emission are not yet 
clear, and detailed contemporaneous multiwave- 
length observations will likely be required to re- 
solve the situation. 

Previous X-ray observations of this source have 
been made with ROSAT, ASCA, BeppoSAX, 
RXTE, XMM, INTEGRAL and Chandra (see [3] 
for a summary), revealing complex behaviour with 
significant variation of flux and photon spectral in- 
dex on both short (~ 1000 s) and longer (monthly) 
timescales. The detection of a TeV signal from LS 
I +61 303 by the VERITAS telescope array in faU 
2006 led us to propose target of opportunity ob- 
servations with the Swift satellite. The resulting 
dataset, combined with some earlier Swift obser- 
vations, provides the most well-sampled long-term 
X-ray /optical/TeV monitoring which has been per- 
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formed on this source. In this paper we focus 
on some initial results obtained with the Swift X- 
ray Telescope. A complementary paper discusses 
the relationship between these observations and the 
TeV emission [12]. We also note that an analy- 
sis reaching similar results and conclusions has re- 
cently been presented by Esposito et al [5]. 

Observations and Data Analysis 

The Swift satellite carries three instruments: the 
Burst Alert Telescope (BAT) operates between 15- 
150 keV and is used primarily as a gamma-ray 
burst trigger; the X-ray Telescope (XRT) is an 
imaging instrument operating from 0.2-10 keV 
with a single photon point spread function of 
18 arcsec (half -power diameter) and a spectral res- 
olution of 140 eV at 5.9 keV (at launch); UVOT, 
the UV/Optical telescope carries optical filters to 
cover the range 170-650 nm. The primary mission 
of the satellite is to investigate gamma-ray bursts 
and their afterglows; however, between gamma-ray 
burst observations Swift is an excellent tool which 
can provide observations of other astrophysical X- 
ray sources, particularly when fast response and/or 
frequent slews are required. 

The XRT observations were made in photon count- 
ing mode on 24 separate days with exposures 
ranging between 700s and 4700s, covering a pe- 
riod from early September 2006 to late Decem- 
ber 2006 and corresponding to a total exposure of 
56 ksecs. The period following the ToO trigger 
(from November 21st to December 24th) is par- 
ticularly well-sampled, with observations approxi- 
mately once every two days. 

The XRT data were processed using the most re- 
cently available standard Swift tools: xrt software 
version 2.0.1, ftools version 6.2, and XSPEC ver- 
sion 12.3.1. The ancillary response files were gen- 
erated with the xrtpipeline task xrtmkarf. The 
spectral analysis presented in this paper used XRT 
data extracted from the 0.3-10 keV energy band, 
thus avoiding instrument response uncertainties 
present below 0.3 keV. 



Results and Discussion 

Figure 1 shows the lightcurve and hardness ra- 
tio as a function of time for the XRT obser- 
vations. The XRT count rate to flux conver- 
sion factor used, obtained using the mean spec- 
tral fit parameters for these data, was 5.7 x 
10~ n crg cm~ 2 s~ 1 /counts s _1 . The systematic 
error associated with any spectral variations is not 
included in the error bars. Figure 2 shows the 
lightcurve and hardness ratio as a function of or- 
bital phase for these same observations, using the 
orbital elements from [2]. Points indicated by the 
same symbol and colour correspond to measure- 
ments taken during the same orbital cycle. Note 
that the lines between points are drawn only to 
help illustrate this - the short timescale variabil- 
ity is such that this simple interpolation can not be 
used to estimate the emission level between obser- 
vations. 

The lightcurve shows clear variability in the flux 
of a factor of 4 - 5. While there are elevated states 
around orbital phases 0.6 - 0.8 and 0.0 - 0.1, it 
is evident that the emission level does not follow 
precisely the same pattern each orbital cycle. The 
variation in the hardness ratio is not statistically 
significant: x 2 /d.o.f. =45/36 for a constant fit, cor- 
responding to a chance probability of 13%. 

The source spectrum for the complete dataset is 
shown in Figure 3. It is well fit by a simple 
power law plus photoelectric absorption model 
(x 2 /d.o.f.=377/370), with no indication of spec- 
tral lines or a thermal black-body component. 
The differential flux at 1 keV is 1.334±g;g^ x 
10~ n crg cm" 2 s~\ the photon index is 1.73±2;ol 
and the neutral Hydrogen column density is 
0.576±g;gi| x 10 22 cm- 2 . Further analysis will ex- 
amine spectral variations as a function of time and 
phase, and explore the possibility of correlations 
between the spectral parameters. 

While there is no doubt that LS I +61 303 is a 
strongly variable X-ray source, previous evidence 
for a consistent orbital modulation of the X-ray 
flux is unconvincing, largely due to the sparseness 
of the data. Our observations reaffirm that the sys- 
tem is a strongly variable X-ray source. Prelim- 
inary analysis reveals no clear evidence for vari- 
ability in the hardness ratio or correlation between 
the hardness ratio and flux. 



30th International Cosmic Ray Conference 



~ 301 1 

I 25r i 



e 20: 



i ♦ 



53980 54000 54020 54040 54060 54080 54100 

Time(MJD) 




53980 54000 54020 54040 54060 54080 54100 

Time(MJD) 



Figure 1: Lightcurve for Swift observations of LS 1 +61 303. The left panel shows the XRT count rate 
(0.2-10 keV) as a function of modified Julian date, with an average flux of 1.4 x 10~ n erg cm~ 2 s . The 
right panel shows the ratio of the 1-10 keV flux over the 0.2-1 keV flux. 
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Figure 2: Lightcurve for Swift observations of LS I +61 303. The left panel shows the XRT count rate 
(0.2-10keV) as a function of orbital phase, the right panel shows the ratio of the 1-10 keV flux over the 
0.2-1 keV flux, also as a function of orbital phase 
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Figure 3: The X-ray spectrum for Swift observations of LS I +61 303. The fit is for a pure power-law with 
photoelectric absorption. 



Chernyakova et al. [3] have claimed evidence for a 
systematic variation of the X-ray spectrum across 
the orbit, with higher X-ray flux producing a harder 
spectrum, as well as for significant variation in the 
column density, Nh- They explain this variation as 
due to the erratic acceleration and cooling of elec- 
trons in the shocks formed between the pulsar wind 
and an inhomogeneous "clumpy wind" from the 
Be star. Detailed, time resolved, spectral analysis 
of our dataset will provide additional constraints to 
this, and other, models for non-thermal emission 
from LSI +61 303. 
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